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ABSTRACT 

A remarkable coincidence has led to the discovery of a family of packings of m? + m — 2 
m/2-dimensional subspaces of m-dimensional space, whenever m is a power of 2. These pack- 
ings meet the "orthoplex bound" and are therefore optimal. 
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1. Introduction 



Let G{m, n) denote the Grassmannian space of all n-dimensional subspaces of real Eu- 
clidean m-dimensional space R™. The principal angles 9i, . . . ,6n G [0, 7r/2] between two sub- 
spaces P, Q £ G{m, n) are defined by 



for i = L . . . , n, subject to u ■ u = v ■ v = 1^ u ■ Uj = v ■ Vj = Q {1 < j < i — 1). We define the 
distancqj between P and Q to be 



In we discussed the problem of finding good packings in G{m,n), that is, for given = 
1, 2, . . ., of choosing Pi, ... , P/v G G{m, n) so that mind(Pi, Pj) is maximized. It was shown 
that for N > m(m -|- l)/2 the highest achievable distance, djs[{m,n), satisfies 



A necessary condition for equality to hold in is that < (m — l)(m -|- 2). An especially 
interesting case occurs when m is even, n = m/2, and A^ = (m — l)(m -|- 2), where we found 
optimal packings for m = 2,4 and 8; that is, packings of 4 lines in M^, 18 2-spaces in R'^ 
and 70 4-spaces in R^. The first is the familiar configuration seen on the British flag (the 
Union Jack), the second is the "double- nine" , a classic configuration from nineteenth-century 
geometry (see the references in |llj and also below), but the third was discovered only after 
a very considerable computer-assisted search. At the time was written we believed that 
there would be no further examples in this series. 

It came as a considerable surprise therefore when we discovered that such packings exist 
whenever m is a power of 2. 

These packings were discovered by a remarkable coincidence. One of us (P.W.S.) had 
discovered a family of groups in connection with quantum coding theory and asked the 
other (N.J.A.S.) for the best way to determine their orders. N.J.A.S. explained to P.W.S. 
that the Magma computer system jH], 0) [H] was ideal for this, and gave as an example the 

symmetry group of above-mentioned set of 70 4-spaces in M®, an 8-dimensional group of order 

^It is shown in that this is a metric, and in fact is essentially the L2 distance between the matrices that 
describe the orthogonal projections onto P and Q. 



cos 9i = max max u ■ v = Ui ■ Vi 

u£P veQ 




2^8! = 5160960. To our astonishment, the first of his groups that P.W.S. tested turned out to 
be (almost) exactly the same group. 

The version of the group that arises from quantum coding in fact has the coordinates in a 
slightly nicer order, and produces the 70 planes as the orbit of the plane spanned by the first 
four coordinate vectors. With the help of our colleague R. H. Hardin we verified that the next 
three groups in the series produced packings meeting the bound in 16, 32 and 64 dimensions. 
Further investigation then produced the general construction given in Section 3. The groups 
are described in Section 2. 

2. The group 

The group Qi that arises from quantum coding theory is a subgroup of the real orthogonal 
group 0{V, M), where V denotes M"*, m = 2*, i > 1, with coordinates indexed by binary i-tuples 
X = (xi, . . . ,Xi) € F ', and F is the field of order 2. Qi is generated by the following 2' x 2* 
orthogonal matrices: 

(i) all permutation matrices tta^ corresponding to affine transformations x i— > Ax + b of F*, 
where A is any invertible i x i matrix over F and 6 € F', and 

(ii) the matrix H = diag{H2, H2, . . . , H2}, where H2 = ;^(]}]!^) (and + denotes +1, — 
denotes —1). 

By multiplying these generators it is easy to see that, for i > 2, Qi contains the matrix 
H' = diag{if4, H4, . . . , H4}, where 



Let Hi be the group generated by the permutations and H'. Then Qi = Hi\JHHi. 

The packings described in Section 3 can be obtained by writing the coordinates in the 
natural lexicographic order and taking the orbit under Qi of the subspace spanned by the first 
2*~^ coordinate vectors (i.e., those in which xi = 0). However, the construction now given in 
Section 3 is a recursive one that no longer explicitly mentions the group. The group is only 
needed in the analysis, where we make use of the fact that it acts transitively on the subspaces. 
In the rest of this section we shall therefore give only a brief discussion of these groups, in 
order to show their connection with the Barnes- Wall lattices. 
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It turns out that Tii and Qi are well-known groups. Tii is the Clifford group CT+(2*) studied 
in 1^, [Sj, jl4j . which in recent years has been used in the classification of finite simple groups 
(see the references in |2j). 7ii \s relevant for the present work because of its connection with 
the Barnes- Wall lattices. 

Although the original paper of Barnes and Wall [21 describes a family of lattices in each 
dimension m = 2* (i > 1), the most interesting lattices are the pair with the highest number of 
minimal vectors (this number is given by the formula displayed in (0J)). We denote this pair of 
2*-dimensional lattices by BWi and BW^. A construction of these lattices using Reed-Muller 
codes is given in [2] and in J2j, p. 234, example (f) (see also [T3]). 

BWi and BW^ are geometrically similar lattices, differing only by a rotation and change 
of scale. When i = 1, for example, we can take BWi to be the square lattice (Fig. 1, solid 
circles), and BW[ to be its sublattice of index 2 (Fig. 1, double circles). In this case the matrix 
D = ^/2H2 acts as an endomorphism sending BWi to BW[. In exactly the same way, the 
matrix y/2H sends BWi to BW^, a geometrically similar sublattice of index 2™/^ (cf. jl2j . 
pp. 240-241). Applying V2H twice sends BWi to 2. BWi. 

(D • (D • (D 
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Figure 1: The first pair of Barnes- Wall lattices, BWi (small circles) and BW[ (large circles). 

Wall showed that for i ^ 3, Tii is the full automorphism group of the lattices BWi and 
BW-. (The case i = 3 is special, since BW3 and BW^ are copies of the root lattice Eg.) The 
group Tii has a normal subgroup E which is an extra-special 2-group of order 2^+^% and TCi/E 
is isomorphic to the orthogonal group Oj(2) = Di{2). The order of Tii is 

i-l 

22i+i . 2i(i-i) (2* - 1) [] (4^' - 1) . (2) 
3=1 

By adjoining the irrational matrix H we obtain the full group ^j, twice the size of 7ii. The 

3 



group Gi also appears in an apparently totally different context in |9| (see the group L defined 
in Eq. (2.13)). 

The way the group Qi arises in quantum coding theory is as follows. The quantum state 
space of i 2-state quantum systems is the complex space C™, m = 2*. Quantum computation 
involves making unitary transformations in this space (see ^01; H)- Some transformations 
may be much easier to realize than others, and it is therefore important to know which sets of 
transformations are sufficient for quantum computation, that is, which sets generate a group 
dense in SU{2^). An interesting set of transformations which generate a finite group are 
the linear Boolean functions on quantum bits (the permutation matrices in our group Qi)^ 
and certain rotations of quantum bits by 7r/2. To obtain the corresponding subgroup of the 
orthogonal group 50(2*), only one rotation is required, which can be taken to be the matrix 
H. 

3. The construction 

We specify a subspace P G G{m, n) by giving a generator matrix, that is, an n x m matrix 
whose rows span P. We will use the same symbol for the subspace and the generator matrix, 
and P-*" will denote the subspace orthogonal to P (or a generator matrix thereof). / denotes 
an identity matrix. 

The construction is recursive. We define a set Qi containing 2^*^^ monomial matrices of 
size 2^-1 X 2^-1 by Qi = {(+),(-)}, 

2^ = {(o +)^^' (o -)®^' (+ o)^^' (- o) ^'-4 ' 

for i > 2. Then Ci is defined by 

Ci = {(+0),(0+), (++),(+-)} , 

Ci = {(/0),(0/),(^^ p)'(o p±)'W); ^eCi_i, QeQij , 
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for i > 2. For example, C2 consists of the 18 matrices 



+ 000 
0+00 

+000 
0+0 



+ + 
+ + 



00+0 
000 + 

+00 
000+ 



+ + 
+ 0- 



0+00 
000+ 

++0 
0+- 



+ 0-0 
+ + 

+ 00 + 
0+-0 



0+00 
00+0 

+-00 
0+- 



+ 0-0 
+ 0- 

+ 0- 
0++0 



++0 
0++ 

+-0 
0++ 



+ 00 + 
0++0 

+ 0- 
0+-0 



(3) 



(The third set of matrices in are the matrices {IQ).) These are generator matrices for 18 
2-spaces in M^. 



Theorem. Let m = 2*, i > 1. The generator matrices Ci define a set of (m — l)(m + 2) = 
22t _|_ 2« — 2 /2m- dimensional subspaces of R™'. The distance between any two subspaces is 



either yjmj^ or yJm/2. 



Proof. The number of subspaces is, by induction, 

2 + 2(2^*"^ + 2*"^ - 2) + 2^*"^ 



22i f 



as claimed. 

Since the recursive definition of the Cj mentions the matrices (/O) and (0/), the coordinate 
positions of C, can be labeled from left to right with binary i-tuples in the natural order, and 
the group Qi then acts by multiplication on the right. It is now easy to find matrices in Qi that 
permute the subspaces transitively. We leave the details to the reader. Therefore, to determine 
the distances between the planes, we may assume that one of the planes has generator matrix 



A 



1 







1 







1 







1 
We recall (cf. |11') that if a second plane has generator matrix 



ci 



C2 



S2 



B 



C3 



S3 
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where + = • • • = + = 1, n = 2*~^, then the principal angles between A and B are 
arccos ci , arccos C2 , . . . , arccos c„ . 

The principal angles between A and (0/) are tt/2 {n times). Between A and the subspaces 

P 0\ fP 

p) °^ U P^ 

they are (n/2 times), 7r/2 (n/2 times); and between A and (/Q) they are 7r/4 (n times). The 



distance from A to any other plane is therefore either y/n/2 or ^/n. ■ 
Since the bound is achieved, this is an optimal packing. 

Together with R. H. Hardin, we are also investigating other families of subspaces that can 
be obtained from the same group. If the initial subspace is taken to be that spanned by the 
first coordinate vector, the orbit consists of the minimal vectors of the Barnes- Wall lattice 
BWi, together with their images under the transformation H, giving a total of 

(2 + 2)(22 + 2)---(2' + 2) (4) 

lines, with minimal angle 7r/4. Taking the plane spanned by the first two coordinates as the 
initial plane, we appear to obtain packings in G{m, 2) containing 

^(2^-l)n(2'- + 2) 

r=0 

planes, with minimal distance 1, for m = 2*, i > 1. 

On the other hand, if the initial subspace is that generated by the first m/4 coordinate 
vectors, we appear to obtain packings in G(m,m/4) containing 

^(m - 2)(m - l)(m + 2)(m + 4) 



subspaces, with minimal distance for m = 2*, i > 2. The first member of this sequence 

is the packing of 24 lines in formed from the diameters of a pair of dual 24-cells. 

We hope to discuss these packings (which appear to be a kind of Grassmannian analogue 
of Reed-Muller codes and Barnes- Wall lattices) in a subsequent paper. 
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